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ABSTRACT 

We analyze observations of the microwave sky made with the Python exper¬ 
iment in its hfth year of operation at the Amundsen-Scott South Pole Station 
in Antarctica. After modeling the noise and constructing a map, we extract 
the cosmic signal from the data. We simultaneously estimate the angular power 
spectrum in eight bands ranging from large {£ ~ 40) to small {i ~ 260) angular 
scales, with power detected in the first six bands. There is a significant rise in the 
power spectrum from large to smaller {£ ~ 200) scales, consistent with that ex¬ 
pected from acoustic oscillations in the early Universe. We compare this Python 
V map to a map made from data taken in the third year of Python. Python 
III observations were made at a frequency of 90 GHz and covered a subset of 
the region of the sky covered by Python V observations, which were made at 40 
GHz. Good agreement is obtained both visually (with a hltered version of the 
map) and via a likelihood ratio test. 
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1. INTRODUCTION 

Since the detection of anisotropy in the cosmic microwave background by the COBE 
satellite, many experiments have measured the angular power spectrum at degree and sub¬ 
degree angular scales (e.g., Netterheld et ah 2002; Halverson et al. 2002; Lee et ah 2001; 
Miller et al. 1999). The Python V data set has sufficient sky coverage to probe the smallest 
scales to which COBE was sensitive, while having a small enough beam to detect the rise 
in the angular power spectrum to degree angular scales, providing a link in Uspace between 
COBE and other recent measurements. 

Python V is the latest of the Python experiments at the South Pole. Dragovan et al. 
(1994), Ruhl et al. (1995), and Platt et al. (1997) describe Python I-III and Rocha, et al. 
(1999) derive constraints on cosmological parameters from these data. Kovac et al. (1997) 
describe the Python IV results. 

The Python V experiment, observations, and data reduction are described in Coble et 
al. (1999). In that paper, we analyzed individual modulations of the data. The modulations 
can be thought of as hlters which have little sensitivity to some of the contaminants in the 
time stream. For example, they have no sensitivity to gradients, which should get a large 
contribution from the atmosphere and from the ground shield. The modulation approach 
also provided a rapid means of compressing a large amount of data (19 Gbytes) into a more 
manageable size. Measurements of anisotropy were reported for eight different modulations 
of the sky signal; the results indicated a sharp rise in the power spectrum. 

In this paper we hud the constraints on the power spectrum due to all of the modulations 
simultaneously. We use the modulations as our starting point, rather than the time stream, 
to take advantage of the contaminant hltering and data compression. We extend the analysis 
of Coble et al. (1999) by accounting for the correlations (in both signal and noise) between 
different modulations. From the modulations we find the best-ht map and its associated 
noise covariance. From this map and its associated covariance matrix we estimate the power 
spectrum simultaneously in eight bands. 

In § 2 we briefly review the instrument and the data set. In § 3 we discuss the estimation 
of the noise matrix. In § 4 we describe how to use this matrix to construct a map and a noise 
matrix for the map. This map is used in § 5 to estimate the angular power spectrum in eight 
bands. In § 6 we check the power spectrum derived from the map with the power spectrum 
derived directly from the modulated data. In § 7 we compare the 40 GHz Python V data 
with the 90 GHz Python HI data (Platt et al. 1997), which covered a subset of the region 
of the sky covered by Python V. We hnd good agreement between the two observations in 
the region of overlap, providing a valuable consistency check. This is another indication of 
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a lack of significant foreground contamination (see also our estimates in Ganga et al. 2002). 
We conclude in 5 8. 


2. INSTRUMENT AND DATA 

We begin with a brief review of the Python V instrument and data, emphasizing the 
terminology used to describe the different subsets of data. More detailed descriptions of the 
instrument can be found in Coble (1999) and Alvarez (1996). A more detailed description 
of the Python V data set can be found in Coble (1999). 

The receiver consists of two focal-plane feeds, each with a single 37-45 GHz HEMT 
amplifier. The two focal-plane feeds of the receiver correspond to two beams at the same 
declination separated by 2.80° on the sky. Each of the two feeds is split into two frequency 
channels near 40 GHz, yielding a total of four data channels. The receiver is mounted on a 
0.75 m-diameter off-axis parabolic telescope, which is surrounded by a large 12-panel ground 
shield. The instrument was calibrated using thermal loads for the DC calibration; the overall 
uncertainty in the calibration of the data set is estimated to be (-1-15%, -12%) in AT. The 
combined absolute and relative pointing uncertainty is estimated to be 0.15°, as determined 
by measurements of the Moon and the Carina nebula {a = 10.73*^, S = —59.65°). The Python 
V beam is well approximated by an asymmetric Gaussian of FWHM 0.91 ^)J;[(?deg X 1.021°; 

01 

deg {az X el), as determined from scans of the Carina nebula and the Moon. Given this beam 
size uncertainty of approximately 0.015°, the band power can move roughly by a factor of 
exp(±/(0.425)(0.015)(7r/180)), only a 3% effect at I = 200. 

Python V observations were taken from November 1996 through February 1997. Two 
regions of sky were observed: the Python V main held, a 7.5° x 67.7° region of sky centered 
at a = 23.18*^, S = —48.58° (J2000) which includes helds measured during the previous 
four seasons of Python observations and a 3.0° x 30.0° region of sky centered at a = 3.00*^, 
S = —62.01° (J2000), which encompasses the region observed with the ACME telescope 
(Gundersen et al. 1995). The total sky coverage for the Python V regions is 598 deg^. 

Both Python V regions are observed with a grid spacing of 0.92° in elevation and 2.5° 
in right ascension, in 345 effective fields. There are 309 unique held positions, but some 
positions are observed at diherent times of the observing season and are thus counted as 
diherent helds for analysis purposes. The telescope is positioned on one of the helds and the 
chopper smoothly scans the beams 17° in azimuth in a nearly triangular wave pattern at 
5.1 Hz. One cycle corresponds to all of the data taken in one back-and-forth scan along the 
sky. A cycle consists of 128 samples along the sky in the given held. A stare is 164 cycles. 
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again centered on the same spot on the sky. One data file consists of roughly ten stares, at 
adjacent helds on the sky. Typically, a hie corresponds to data taken over hve to ten minutes 
(depending on how many stares it contains). The telescope remains on this set of helds for 
roughly 13 hours, so any set of helds is typically observed in about one hundred consecutive 
hies. This terminology is summarized in Table 1 and illustrated in Figure 1. 

In software, the data are modulated such that the spatial responses are cosines apodized 
with a Hann window. In order to take advantage of the large sky coverage of Python V, 
which allows us to probe large angular scales, we also use an additional cosine modulation 
which was not apodized by a Hann window. Data taken during the right and left-going 
portions of the chopper cycle are modulated separately, to allow for cross-checks of the data. 
Sine modulations are not used in the analysis because they are anti-symmetric and are thus 
sensitive to gradients. The modulated data in a given stare is a linear combination of the 
samples: 

128 

Dima 'y ^ ^msdisa- ( 1 ) 

5 = 1 

The index i = 1,...,690 labels the held and feed; m = 1,...,8 labels the modulation; 
a = 1,..., Nf labels the hie which looks at a given set of helds; s = 1,..., 128 indexes the 
sample number; d is the unmodulated data which has been co-added over all cycles in a 
stare. 

A chopper synchronous ohset, due to dihering amounts of spillover, is removed from 
each data hie by subtracting the average of all stares in a hie. This is not just a DC ohset; 
there is an ohset removed for each modulation. The chopper synchronous ohset is discussed 
in detail in Coble (1999), but typical values are 100 to 200 /iK for each modulation, stable 
over a timescale of ~ 5 hies. 

When the data are binned in terrestrial azimuth, a periodic signal due to the 12 panels 
of the ground shield is evident, especially in the lower-£ modulations. This signal of period 
30° is ht for an amplitude and is subtracted. Removal of the ground shield ohset has less 
than 4% ehect on the hnal angular power spectrum because when the data are binned in 
RA, the ehect averages out. The ground shield ohset is discussed in detail in Coble (1999), 
but typical values for the signal amplitude are less than 100 /iK. 

Both the chopper synchronous ohset and the ground shield ohset subtractions are ac¬ 
counted for by adding a constraint matrix, C^, to the noise matrix (Bond, Jahe, & Knox 
1998). The precise form of these matrices is given in Coble (1999). Their impact on the hnal 
result is minimal since they serve to remove only a handful of modes from the analysis. In 
the Coble et ah (1999) analysis, the chopper synchronous ohset and the ground shield signal 
were removed from the data, but the ground shield constraint matrix was not included in 



5 


the analysis. The constraint matrix for the chopper synchronous offset was included in that 
analysis. 

After the data have been modulated and offsets removed, the right- and left-going data, 
which have been properly phased, are co-added, as are data from channels which observe 
the same points on the sky. Data pointing at a held i from all hies are averaged to form 


Dim / Dima- 

iV f 

J n = ^ 


This hnal data vector has 5520 (= 345 helds x 2 feeds x 8 modulations) components. The 
next section describes our modeling of the noise properties of these data. 


3. NOISE MODEL 

Accurate modeling of the noise is often one of the most difficult tasks in CMB analysis. 
The noise model we develop below enables us to estimate the angular power spectrum in 
eight bands simultaneously. See Coble (1999) for a more detailed discussion of the noise 
modeling of the Python V data. In Coble et ah (1999), we modeled the noise only for 
individual modulations. The noise model described here also models the cross-modulation 
terms, allowing us to include cross-modulation correlations in the power spectrum analysis. 
The noise level for the Python V data is < ImKs^'^^. 

Our noise model assumes that the covariance between helds taken with diherent sets 
of hies is negligible (as in Coble et ah 1999) because of the chopper ohset removal and 
because of the long time between measurements. An analysis comparing the noise estimated 
on diherent timescales indicates that Python V noise is dominated by detector noise and is 
Gaussian. 

Since many diherent hies look at the same held on the sky, there is a simple way to 
estimate the noise covariance matrix. We hrst estimate the noise matrix via: 

where again i,j index the diherent helds, m,m' the eight modulations, and we sum over all 
Nf hies which observe the helds of interest. However, since there are typically only 100 hies 
for each held, the sample variance on the noise estimate is 1/(100)^/^ or 10%, which will 
severely bias estimates of band power. Hence we do not use this naive estimator. 






Dim){D 


jm'a Djm' 
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To obtain a better estimate of the noise, in Coble et al. (1999), we averaged the variances 
for each set of hies and then scaled the off-diagonal elements of the covariance to the average 
variance in a given set based on a model derived from the entire Python V data set. In 
that paper, was computed for each individual modulation, i.e., with m = m! only. 

Several consistency checks were performed showing that the hnal noise model for the single 
modulation analysis was a good one. 

In this cross-modulation analysis, we initially extended the method in Coble et al. 
(1999) to account for cross-modulation terms in i.e., terms with m ^ m!. To test 

this noise model, we constructed = D*(C'^ -|- C‘")“^D, for each observing set (which 
typically includes of order ten helds observed ~ 100 times each). There is very little CMB 
signal in any one set, so we expect y^/dof to be close to one. The results fail this test, 
indicating that a better model of the cross-modulation noise is necessary. 


To go beyond the initial estimators, we assume the cross-modulation noise matrix factors 


as 


^ ij mm' ^ mm' ^ ij ’ 


( 4 ) 


where describes the cross-modulation correlations and the held-field correlations. 
The cross-modulation correlations are derived from the sample-space covariance matrix C®: 


^ mm' 


128 


MmsCggiMm's'- 

s^s'=\ 


( 5 ) 


The matrix C® describes the noise in the timestream as a function of chopper sample s. To 
clarify, is a 128 x 128 matrix, is a 8 x 8 matrix, and is a 690 x 690 matrix for 
Python V. Models for and are needed in order to construct C^. 

If we assume that C® depends only on chopper sample separation As = s — s', it can 
be computed from the following function: 

/(As) = (6) 

Ns^ 

where Ns is the number of samples and dg is the unmodulated data. For example, the Cf 2 
component is given by /(As = 1). In order to compute /(As), a chopper synchronous offset 
is first subtracted from the raw data. Then /(As) is calculated for each channel, cycle and 
stare in a hie. /(As) is then averaged over cycles, stares, and hies. Figure 2 shows /(As) 
for each channel in one of the sets. 


With this model for the sample correlation function, it is now straightforward to compute 
for each set and channel following equation (5). matrices for channels which look at 
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the same point on the sky and for right and left-going chopper data are averaged, yielding 

matrices for both feeds in each set. As an example, for one set and feed is shown 
in Figure 3. 

In order to get a simple form for C^, the held correlation matrix, we ignore the cor¬ 
relations between the two feeds and assume the correlation between helds i and j come 
only from the chopper offset subtraction. We investigated several similar noise models and 
found that these assumptions do not change the single modulation angular power spectrum 
signihcantly, so we assume is of this form for the cross-modulation analysis. 

Finally, since this noise model is derived from sample to sample huctuations, it is larger 
than the corresponding noise derived from the co-added data by a factor of ~ 10^, so 
must be normalized to the variance in the co-added data for each set. Since accounts 
for the relative normalization of all of the modulations, must be normalized to the 

variance in only one modulation of the co-added data for each set. We normalize to modu¬ 
lation 8 because we expect the higher order modulations to be least affected by the ground 
shield. Figure 4 shows the y^/dof for each set using the hnal cross-modulation noise model, 
indicating a good final noise model. 

As another check on the noise matrix used in the cross-modulation analysis, single 
modulation band powers were computed using the components of These are 

consistent with the band powers given in Coble et al. (1999). 


4. MAPS 

We want to estimate the power spectrum from the Python V data. Ideally, given the 
non-circular beam, this should be done directly from the modulated data. This would require 
us to form the likelihood function from the covariance matrix for the N = 5520 data points. 
Inversion or Cholesky decomposition of matrices are processes so computational demands 
are signihcantly alleviated by creating a map (with N = 1666 at highest resolution) from the 
modulated data and then estimating the power spectrum from the map. This technique was 
used in an analysis of the MSAM-I experiment (Wilson et al. 2000), for which the power 
spectrum estimated using the map is consistent with the power spectrum estimated directly 
from the modulated data. The inversion problem of map making and the circular beam 
assumption used in it does call for cross-checks and verihcation against likelihood analyses 
of the modulated data. The results of such tests are summarized in § 6. 

The data can be expressed as: 


D = MT -h n 


( 7 ) 



with noise covariance matrix {nn) = N = + C*". As mentioned above, the data vector D 

has 5520 elements. The matrix M describes the experimental processing of the underlying 
temperature held; it is equal to the modulations with an index corresponding to each pixel 
at which we estimate the temperature T. Given the modeling of the data as in equation (7), 
the minimum variance estimator for T is 

f = NMN-^D. (8) 

This estimator will be distributed around the true temperature due to noise, where N, the 
noise covariance matrix for the map, is given by 

N= ((f-T)(f-T)) = (9) 

The inversion in equation (9) is singular, so it is performed via Singular Value Decom¬ 
position. It is obvious which modes are singular and should be neglected. We have tested 
various thresholds and found no change in the results. The pixels in the map are 30' in RA 
which corresponds to about 20' on the sky. Coarser grids gave similar results for the band 
powers; as we will see, there is little sensitivity to modes with ^ > 200, so 20' (a third of the 
beam size) is more than adequate. 

Another advantage of the map basis is that the theory covariance matrix is simple to 
compute. In the map basis, the theory covariance matrix simplihes to 

Cl =< T,T, >= ^^P,(cos%)e-^'"'C',, (10) 

i 

where i and j now refer to map pixels. Pi is a Legendre polynomial, and Qij is the angular 
separation between points. We take a = (o'azO'ei)^'^^ = 0.425 x 0.96°. Taking the beam to 
be circular will not change the band powers signihcantly (see § 6). From equation (10), the 
window functions in the map basis only depend on the angular separation Qij and not on any 
of the details of the observing strategy. This is a smaller basis than that used in the analysis 
of the modulated data which accounts for beam non-circularity, described in § 6. Indeed, 
one way to think of a map is that it is the linear combination of the data for which the signal 
(and therefore its covariance) is nearly independent of the specihc experimental observing 
strategy. The noise covariance for the map (eq. [9]) accounts for all of the experimental 
processing and the constraints. 

Although we are primarily interested in the map as a vehicle on the road to the power 
spectrum, it can also be Wiener-hltered to produce a realistic image of the sky. Wiener- 
hltered maps of both of the Python V regions are shown in Figure 5. We use the unhltered 
map for power spectrum estimation. The map serves another useful function apart from its 
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use for the power spectrum. One can use maps to compare different data sets that were 
processed in completely different manners. In § 7 we present a visual comparison of Python 
III and Python V. First, though, let us compute the power spectrum. 


5. ANGULAR POWER SPECTRUM FROM THE MAP 


We now use the map to estimate the CMB anisotropy power spectrum. Because the 
observations are far short of full sky coverage, we cannot determine individual C^’s. Instead, 
we parameterize the theory covariance matrix, C'^, with the power spectrum, C = ^(^ + 
l)Q/(27r), broken into bands of Q, denoted by a 


= = | ^ ( 11 ) 

a ^ 

so that within each band a, Ci = Ca = constant and iminio) and imaxici) delimit the range 
of band a. We use eight contiguous bands of equal width, as given in Table 2. 

Since the CMB anisotropy appears to be Gaussian on the angular scales probed by the 
Python V experiment (Park et al. 2001; Wu et al. 2001; Shandarin et al. 2002), we can 
in principle use the theory covariance matrix for the map (eq. [10]) together with the map 
noise matrix, N, and the pixelized map data, T, to form the full likelihood function: 

C = (27r)“^'^^det(C)~^'^^exp(—x^/2) (12) 


where = T*C~^T and C = + N. We can then hnd the Ca which maximize it 

by conducting a direct, grid-based search in the full eight-dimensional parameter space. 
In practice, this is of course unfeasible because it would require of order ten likelihood 
evaluations in every dimension of parameter space. The likelihood function computation 
requires an inversion and a determinant of a large matrix (in our hnest pixelization, 1666 x 
1666), so it is certainly impractical to attempt this 10® times. 

Instead, we use the quadratic estimator (Bond, Jaffe, & Knox 1998; Tegmark 1997) to 
hnd the maximum likelihood band powers and their errors. Dehning 


A. 


a,ij 



dCa 



(13) 


where C is the full theory plus noise covariance matrix, the Fisher Matrix which describes 
the errors is 


Fab = 



OCt 



(14) 
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and the quadratic estimator is 

c„ = cr + If-/(ta^t - Tr c-‘^ ), (15) 

We start from a flat spectrum (e.g. all = 1000 and iterate four times. Convergence 
to well within the size of the error bars is usually reached by the second iteration. 

The band temperature AT£{= results of the likelihood analysis are shown in Figure 
6 and given in Table 2. The AT^ error bars here account only for the statistical uncertainties, 
and in particular, do not account for the calibration or beam uncertainties. The Fisher matrix 
is given in Table 3. We emphasize that the values of the angular power spectrum differ from 
those in Coble et ah (1999) because we are including more information in this analysis: 
the cross-modulation correlations. Again, single modulation band powers computed using 
the components of are consistent with the band powers given in Coble et al. 

(1999). 

6. Comparison of Band Temperatures from Modulated Data and Map 

The map has a much smaller basis {N = 1666 at the hnest pixelization) than the 
modulated data {N = 5520) and hence allows for speedier likelihood analysis. However, 
map making is an extra step that needs verihcation. In particular, the map analysis has to 
assume a circular beam when constructing the theory covariance matrix C'^. 

In the modulated data basis, the beam corresponds unambiguously to the measured 
beam response function. The non-circular (elliptical) Python V beam is an additional com¬ 
plication for C'^ computations. Souradeep & Ratra (2001) develop computationally rapid 
methods for computing for experiments with non-circular beams. The constant eleva¬ 
tion scans of the Python V experiment allow us to exactly incorporate the effects of beam 
non-circularity, without recourse to any approximation. 

The larger size of the modulated data basis makes the 8 band likelihood analysis de¬ 
scribed in § 5 computationally expensive. We therefore choose to compare the map basis 
and modulated data basis results using a simplihed analysis that accounts for the cross¬ 
correlation between modulations in a limited manner. This likelihood analysis estimates 
the band temperatures in each of the 8 Aspace bins while holding hxed the other 7 band 
temperatures at the central values obtained in § 5. 

Table 4 compares the band temperature estimates in the map basis and the modulated 
data basis. The two sets of results agree to 0.5 a. The differences become larger for higher 
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^ bins and possibly arise from non-circular beam effects. Souradeep & Ratra (2001) show 
that non-circular beam effects become more important above the ^ value corresponding to 
the inverse beamwidth. 

The Python V band powers can be used in combination with the results of other ex¬ 
periments to test for consistency and constrain cosmological parameters. This can be done 
in a way that accounts for the non-Gaussianity of the band power uncertainty by using the 
offset log-normal form for the likelihood given in Bond, Jaffe, & Knox (2000): 

-2\nC = Y,(Zl-Zt)Mi(Zl-Zt) (16) 

a,b 

where 

Z^ = ln{Ca + Xa), (17) 

Z*^ln(5^G/iV, + a;J, (18) 

\ l£a / 

and 

= Fab {Ca + Xa) (Cfe + Xb) ■ (19) 

Again a, b denote bands. We have approximated the band-power window function as a 
tophat with width Ni = 30. 

We have £t the Xa and Ca parameters of this form to our one-dimensional likelihood 
curves, as directly evaluated from the modulated data in Table 4. The Fisher matrix comes 
from the quadratic estimator applied to the maps. Table 5 gives the parameters of the offset 
log-normal analytic fits to the band power likelihoods. 

Table 6 compares band temperatures estimated with and without the circular beam 
approximation, from single modulation analyses where the correlations between modulations 
are ignored (as in Coble et al. 1999). The last column corrects the results obtained by Coble 
et al. (1999) for a systematic underestimation of the error bars by a factor of log^Q e. These 
results also use the non-circular beam and do not use the flat-sky approximation (Souradeep 
& Ratra 2001). The effect of the circular beam approximation on the Python V power 
spectrum is minimal, but would be greater for an experiment with higher S/N. 

7. COMPARISON WITH PYTHON HI 

Using the technique of § 4, maps of the sky are constructed from the Python III data. 
We compared the maps in two different ways. First, we decomposed each map into its signal 
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to noise eigenmodes. Keeping all the modes results in little useful visual information since 
most of the features in such a map are noise. Therefore, we excluded all modes with S/N 
less than 1.7; stopping at S/N = 1 retains too much noise. The resulting maps are shown 
in Figure 7. Since Python III has higher S/N than Python V, it retains many more modes. 
Therefore, not all the features seen in the Python III map should be visible in the Python 
V map. However, structures found in the Python V map are evident in the Python III map, 
implying that Python III and Python V are consistent with each other. 


While the visual comparison is quite useful, it is difficult to judge the signihcance of the 
agreement in this manner. For a more quantitative comparison we use the (3 test of Knox 
et ah (1998). This statistic has a number of possible interpretations, one of which is that 
it is the log of the “probability enhancement factor”. That is, it tells us how much more 
probable the data sets 1 and 2 are viewed jointly, as opposed to disjointly: 




In 


P(Ti,T2) 

P(Ti)P(T2)- 


( 20 ) 


This can be re-written in terms of the likelihood function: 


P = ln£(Ti,T 2 ) - ln£(Ti) - ln£(T 2 ) (21) 

The joint likelihood for the two data sets uses the likelihood equation with the data and 
the noise covariance being a concatenation of those from the two data sets. It also uses the 
theory covariance between the two sets. 

We hud (3 = 6.9, which means that the data is e®'® = 992 times more probable viewed 
jointly than disjointly. This shows that the data sets have significantly more in common 
than they would if they were unrelated to each other. 

We can also examine how likely this value of (3 is, under these two different assumptions. 
The hrst assumption is that the data sets are related to each other exactly as we expect due 
to their locations on the sky, and our inference of the signal power spectrum from the Python 
V data. With assumption 1 we hud that {{3) = 11.7 ±4.3. If instead we assume that the two 
data sets are completely unrelated (perhaps because one is completely contaminated), then 
(/5) = —19.9 ±9.4. We see that f3 differs by 1.1 standard deviations from the expected value 
and 2.8 standard deviations from the value expected in the absence of eross-correlations. 

The f3 statistic is model-dependent, but we found that it changed by less than 15% as 
we varied the amplitude of the assumed power spectrum by amounts consistent with the 
error bars and as we adjusted the calibration by ±30%. 
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8. CONCLUSIONS 

The Python V experiment densely samples 598 square degrees of the microwave sky and 
constrains the CMB anisotropy angular power spectrum from £ ~ 40 to £ ~ 260, showing 
that power is increasing from large to smaller {i ~ 200) angular scales. The noise matrix 
constructed in § 3 enables us to simultaneously estimate the angular power spectrum in eight 
bands. The power spectra estimated from the map and directly from the modulated data 
are consistent. The rise seen in Figure 6 is characteristic of acoustic oscillations in the early 
Universe. A number of other measurements also indicate such a rise in power (e.g., Ganga, 
Ratra, & Sugiyama 1996; Netterheld et ah 1997; de Oliveira-Costa et ah 1998; Torbet et al. 
1999; Podariu et al. 2001; as well as experiments mentioned in § 1). Python V extends to 
larger scales (lower i) than these, to the smallest scales to which COBE was sensitive. 

The Python III and V experiments differ in signihcant ways, including frequency, re¬ 
ceiver, year, and noise properties. Nevertheless, the maps and the (3 test in § 7 indicate that 
they both detect similar signals, a rare and very valuable consistency check and conhrmation. 
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Fig. 1.— Observing sets. Each set is observed for approximately a day (leaving ~ 13 hours 
of good data) before moving on to the next set. Neighboring sets in the main Python V 
region overlap by 1 held. Circles represent the helds and arrows point to the end helds in 
each set. In addition to the sets shown, sets ib, jb, kb, and lb were observed with a scan 
pattern of 5 helds per hie. 
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Fig. 2.— Correlation as a function of sample separation, /(As), for all four data channels 
in one of the sets of hies. The noise is a combination of atmospheric and instrumental noise. 
Channel 5 is our most sensitive channel. Channel 3 is a dark channel. These noise levels 
correspond to < ImKs^'^^. 
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Fig. 3.— for one set and feed. Elements that are more than two modulations apart are 

relatively uncorrelated. 
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Fig. 4.— x^/dof for each set using the final noise matrix. Typical number of dof is 96 (= 
13 stares x 8 modulations less 8 constraints), though some are much smaller. The y^/dof 
close to 1 indicates that the hnal cross-modulation noise model is a good one. 
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Fig. 5.— Wiener-filtered CMB map for the Main and South Python V regions. The maps 
are plotted with the same size and temperature scales. The unfiltered map was used for 
power spectrum estimation. 
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Fig. 6.— Angular power spectrum from the cross-modulation analysis. Band powers are 
simultaneously estimated in 8 t'-space bands. Horizontal bars indicate the width of the 
bands. Errors are la statistical only and do not include calibration or beam uncertainties. 
The last two bins do not show 2a detections. The corresponding 2a upper limits on AT^ are 
174 and 211 pK, respectively. 
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Fig. 7 .— Comparison of Python V and Python III maps. The maps are made using only 
the highest S/N modes. Structures found in the Python V map are evident in the Python 
III map, implying that Python III and Python V are consistent with each other. 
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Term Definition 

Field Center position for an observation. 

Cycle One back-and-fortli scan along the sky centered on one 
field. Consists of 128 samples. 

Stare 164 consecutive cycles, again centered on the same held. 

File Approximately 10 consecutive stares. Each stare is 
centered on a new adjacent held. 

Set Approximately 100 consecutive hies (about 13 hours of 
data taking) of the same ~ 10 helds. 


Table 1; Data subset terminology. See Figure 1 for a visual representation of the sets of 
helds. 




Bin 

i 

at, (/iK) 

1 

44^!| 

221^ 

2 

75t!t 

24^7 

3 

loeili 

3419^ 

4 

i37ili 


5 

168^11 

611}? 

6 

i99iii 

771?° 

7 

230il| 

00 o 
+ 1 

O 

8 

26iili 

nq+71 


Table 2: Python V final angular power spectrum. This comes from the cross-modulation 
analysis. Note- Band powers are simultaneously estimated in 8 £-space bands. Error bars 
are la statistical only and do not account for the calibration or beam uncertainties. The 
last two bins do not show 2a detections. The corresponding 2a upper limits on AT^ are 174 
and 211 /iK, respectively. 
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1 

2 

3 

4 

5 

6 

7 

8 

1 

1.00 

-0.164 

0.012 

-0.008 

-0.003 

-0.003 

-0.002 

-0.002 

2 

-0.164 

1.00 

-0.211 

0.019 

-0.011 

-0.004 

-0.002 

-0.004 

3 

0.012 

-0.211 

1.00 

-0.217 

0.024 

-0.012 

0.000 

-0.014 

4 

-0.008 

0.019 

-0.217 

1.00 

-0.228 

0.030 

0.005 

-0.064 

5 

-0.003 

-0.011 

0.024 

-0.228 

1.00 

-0.236 

0.017 

-0.089 

6 

-0.003 

-0.004 

-0.012 

0.030 

-0.236 

1.00 

-0.361 

-0.003 

7 

-0.002 

-0.002 

0.000 

0.005 

0.017 

-0.361 

1.00 

-0.385 

8 

-0.002 

-0.004 

-0.014 

-0.064 

-0.089 

-0.003 

-0.385 

1.00 


Table 3: Fisher Matrix FJ/{^ 








Bin 


AT, (/iK) AT, (/iK) 

(Map) (Modulated Data) 

1 

4411^ 

25^^ 

25^^ 

2 

7511^ 

23^^ 

22^6 

3 

106^}^ 

34^^ 

32+® 

4 

137^}^ 


50^^ 

5 

168^}^ 



6 

199^}^ 

74^^ 

671^^ 

7 

230^}^ 


q+53 

^-9 

8 

26 u;= 

46^1 

70+39 

<U_7o 


Table 4: Comparison of angular power spectra from the map and from the modulated data. 
Note- These are from analyses that account for the cross-correlations between modulations 
in a more limited manner than the full simultaneous band power estimation of Table 2, as 
described in the text. Error bars are la statistical only. The last two bins do not show 2a 
detections. The corresponding 2a upper limits on AT^ are 112 and 202 juK (from the map) 
and 104 and 181 juK (from the modulated data), respectively. 




Bin 

a;(/iK2) 

Ca(/iK2) 

1 

100 

620 

2 

200 

500 

3 

500 

1050 

4 

5000 

2600 

5 

5000 

3150 

6 

6000 

4570 

7 

20000 

0 

8 

60000 

5000 


Table 5: Parameters of the offset log-normal analytic form for the band power likelihood. 




Modulation £ 

AT, (/iK) 
(circular beam) 

AT, (/iK) 
(elliptical beam) 

1 

50^^^ 

231^ 

231" 

2 

74^5^ 

241^ 

241^ 

3 

1081^? 

30^^ 

30l^ 

4 




5 

^72+tl 



6 

203^38 

102^24 

qci +22 

^'-’-22 

7 

oqq+40 



8 

264^^^ 

oT 

q+78 

U-0 


Table 6: Comparison of single modulation angular power spectra with and without the cir¬ 
cular beam approximation. Note- These are from analyses that ignore correlations between 
the modulations. The last column corrects the results obtained by Coble et al. (1999) for a 
systematic underestimation of the error bars by a factor of log^Q e. Error bars are la statis¬ 
tical only. The last two bins do not show 2a detections. The corresponding 2a upper limits 
on AT^ are 158 and 165 /iK (for the circular beam) and 142 and 143 juK (for the elliptical 
beam), respectively. 




